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I. INTRODUCTION 


The nonlinear Schrodinger (NLS) equation, as one of the universal equations that describe the 
evolution of slowly varying packets of quasi-monochromatic waves in weakly nonlinear dispersive 
media, has been very successful in many applications such as nonlinear optics and water waves^^. 
The NLS equation is integrable and can be solved by the inverse scattering transform^. However, 
in the regime of ultra-short pulses where the width of optical pulse is in the order of femtosecond 
s), the NLS equation becomes less accurate. Description of ultra-short processes requires 
a modification of going beyond the standard slow varying envelope approximation (SVEA) . Re¬ 
cently, Schafer and Wayne derived a short pulse (SP) equation 

Uxt U ~{u ^xx ( 1 ) 

D 

in attempting to describe the propagation of ultra-short optical pulses in nonlinear media:^. Here, 
u = u{x,t) is a real-valued function, representing the magnitude of the electric field, the subscripts 
t and X denote partial differentiation. It has been shown that the SP equation performs better than 
NLS under this case^. 

Apart from the context of nonlinear optics, the SP equation has also been derived as an inte¬ 
grable differential equation associated with pseudospherical surfaces^. The SP equation has been 
shown to be completely integrable^—. The loop soliton solutions as well as smooth soliton so¬ 
lutions of the SP equation were found ii>l^. The connection between the SP equation and the 
sine-Gordon equation through the hodograph transformation was clarified, and then the A^-soliton 
solutions including multi-loop and multi-breather ones were given inl^dl by using Hirota’s bilniear 
methodA^. An integrable discretization of the SP equation was constructed by the authors in^, and 
its geometric interpretation was given in^. 

A major simplification made in the derivation of the short pulse equation is to assume that 
the polarization is preserved during its propagating inside an optical fiber. However, this is not 
always the case in practice. For example, we have to take into account the effects of polarization or 
anisotropy in birefringent fibers^. Therefore, an extension to a two-component version of the short 
pulse equation is needed in order to describe the propagation of ultra-short pulse in birefringent 
fibers. In fact, several integrable coupled short pulse have been proposed in the literature^^— . The 
bi-Hamiltonian structures for several coupled short pulse equations were obtained ii^. 

In the present paper, we are concerned with the integrable semi-discretization of a multi- 
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component short pulse (MCSP) equation 
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^ ^ CjkUjUk I Ui^x I ) 

j<k<n ) / „ 


^ ^ i = 1,2, • • • ,n, (2) 

where the eoeffieients Cj^ are arbitrary eonstants with symmetry c^fc = Cfcj. Eq. (|2l) was proposed 
by Matsuno^^ through Hirota’s bilinear method, meanwhile, multi-soliton solution was given as 
well. 

The remainder of the present paper is organized as follows. In Seetion 2, the MCSP equation is 
briefly reviewed. We provide its iV-soliton solution in an alternative pfaffian form and prove it by 
the pfaffian teehnique. In Seetion 3, by using a Baeklund transformation of the bilinear equations 
and defining a diserete hodograph transformation, we eonstruet a semi-diserete analogue of the 
MCSP equation. Meanwhile, iV-soliton solution in terms of pfaffian is provided and proved. In 
Seetion 4, we investigate in detail the one- and two-soliton solutions to the semi-diserete eomplex 
short pulse equation, whieh ean be reduced from the MCSP equation. The paper is eoncluded by 
several remarks in Seetion 5. 


II. REVIEW OF THE MULTI-COMPONENT SHORT PULSE EQUATION AND ITS 
MULTI-SOLITON SOLUTION 


It was shown by Matsuno in^^ that the SP equation ([T]) is derived from bilinear equations 

■! = \ (f - f 


through transformations 


DsDyf-f'=l(f-P) , 


/' 


(3) 




(4) 


Here D is ealled Hirota H-operator defined by 

f d d \'^ f d d 

o-D-f /(,.. 

Recently, in view of the faet that the SP equation ([T]) ean also be derived from another set of 
bilinear equations 

DsDyf -9 = fg, 

Dsf ■ / = ^^/^ 


(5) 
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through transformations 


X = y - 2{\nf)s, t = s, 


u = 


/’ 


( 6 ) 


Matsuno^^ constructed a multi-component generalization of the short pulse equation O based on 
a multi-component generalization of bilinear equations dS]), which reads 


DsDyf ■ gi = fgi, i = 1,2, 


■ f =2 Cjkgjgk 


(7) 


l<j<k<n 

Remark 1. The set of bilinear equations (12]) is actually obtained from a 2-reduction of the KP-Toda 
hierarchy, which basically delivers only two tau-functions out of a sequence of the tau-functions. 
Furthermore, when these two tau-functions are made complex conjugate to each other, the bilinear 
equations (12|) is converted into the sine-Gordon equation fys = sin f via a transformation f = 
2iln(f'/f), which is further converted into the SP equation ^by a hodograph transformation. 


Remark 2. Ir^, Hirota and one of the authors have shown that both the bilinear equations (12|) 
and (O derive the sine-Gordon equation. Furthermore, the relations between tau-functions, which 
read 

/ = /'/, g = 2iDj'-f, (8) 

were also presented. 


Remark 3. As mentioned previously, Eqs. (12|) originate from 2-reduction of single component 
KP-Toda hierarchy, whereas, Eqs. (l5|) come from (1 -f l)-reduction of two-component KP-Toda 
hierarchy. Since they both belong to of the Euclidean Lie algebrc^, it is natural that they 
both derive the SP equation. However, the latter can be easily extended to (1 -I — ■ -f 1)-reduction 
of multi-component KP-Toda hierarchy, which gives rise to the multi-component generalization of 
the short pulse equation. 


In what follows, we will briefly review how the bilinear equations ([7]) determines a multi- 
component generalization of the SP equationdH). Dividing both sides by p, the bilinear equations 
([7]) can be cast into 


(pj +2|(ln/),„ = |, i = l,2,.--,n, 

n 1 djdk 

l^j<k<n 


( 9 ) 
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Introducing a hodograph transformation 


X = y - 2{\nf)s, t = s, 


( 10 ) 


and a dependent variable transformation 


9i 


we then have 


= (i = ,n), 


— = -2(\nf),s = -- CjkUjUk, — = l- 2{\nf)sy, 


l<j<k<n 


dy 


whieh implies 


dy = p ds = dt- 


vl<j<fc< 

by letting 1 - 2{ln f)sy = p~^. 

Notiee that the first equation in @ ean be rewritten as 


^ ^ CjkUjUk I dx 

j<k<n / 


or 


whieh is eonverted into 


dx dt 


jj =(l-2(ln/),,)|, 


, 9i \ _ 9i 


^ ^ CjkUjUk I dx I 

'j<k<n / / 


Vjj^ Ui , 


( 11 ) 


( 12 ) 


(13) 


vl<j<fc<r 

by the eonversion relation (fT^ . Obviously, Eq. (fTSl) is nothing but the MCSP Q. 

Next, we give an alternative representation of the A^-solution to the MCSP equation Q in the 
form of pfaffians. To this end, let us define a elass of set B^, /i = 1, 2, • • • , n, whieh satisfies the 
following eondition. 


= if /r 7^ V, = {&i, &2, • • • , &27v}- 


Then we define the elements of the pfaffians (others not mentioned below are zeros) 

Pf (a„ ak) = , pf hk) = 6,,k , 

Pj + Pk 

Ff{b,,bk) = 7 _2 , (bj e B^, bk e B,) , 

^Pj -Pk 


(14) 

(15) 
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Fi{di,ak) =pie^'^, 


( 16 ) 




1 bj G 

0 bj ^ Bfj, 


( 17 ) 


Here j,k = 1, 2, • • • , 2N, /r, z/ = 1, 2, • • • , n, = pji/ + pj s + .^*0 and I is an integer. 
By defining the elements of the pfaffians, we ean give the pfaffian solutions satisfying bilinear 
equations (|7]). 

Theorem 1. The bilinear equations ([71) have the following pfaffian solution 


/ — Pf(ai, • • • , a2Nt bi, - ■ ■ , &2Af), 

9% = Pf(<^0! Aj Ol) • ■ ■ 5 tl2N, bl, - ' ' 5 Aiv) 5 


(18) 

(19) 


where z = 1, 2, • • • , n and the elements of the pfaffians are given in Eqs. 47?l)-07l). 


Proof Sinee 


A 

dy 


Pf(aj, af) = {pj - Pfc)e^^'+^'= = Pf(cio, A, Oj, Ofe), 


d 

—Fi{aj,ak) = {p^^ = Pf(d_i, do, Oj, afc), 


—Pf(aj,afc) = {p~^ - = Pf(d_2, do, a^-, Ofc), 


a" 

dyds 


Pf(aj, Ok) = {pjPk - PkPj = Pf(d_i, di, a*, Oj ), 


where Pf(d/, dm) = 0 for integers I and m, we then have 

|^ = Pf(do,di,---), 

dy 


df 

^=Pf(d_i,do, 


ds"^ 

d^f 

dyds 


— Pf(d_2, do, • • •), 


= Pf(d_i,di,---). 


6 





Here Pf((io, ai, • • • , a2N, Wr ■ ■ , ^n) is abbreviated by Pf((io, di, • • ■), so as other similar 
pfaffians. 

Furthermore, it ean be shown 

■ 2N 

5^(-iyPf(do,a,)Pf(A 


dgi _ d 
dy dy 


-'D ) 


-i=i 


2N 


= [{dyPf{do, a,)) Pf (A, •••, a,■,••■) + Pf{do, a,)9,Pf(A, •••,%,•••)] 

i=i 

2N 

= [Pf(di, a,)Pf(A, • • • , di, • • •) + Pf(^^o, a,)Pf(A, do, A, • • • , a„ • • •)] 

i=i 

= Pf((ii, A, • • •) + Pf(do, A,c?o, A, • • •) 

= Pf(A,A,---)- 


Here dj means that the index j is omitted. Similarly, we ean show 


■ 2N 

5^(-iyPf(d_i,aAPf(A,-- - ,%,•••) 

.i=i 

2N 

= [(dyPiid-u ay) Pf(A, • • • , di, • • ■) + Pf(d-1, a,)APf(A, •••,%,•••)] 

i=i 

2N 

= Y.{-iy [Pf(do, aAPf(A, • • • , a„ • • ■) + Pf(d-i, aAPf(A, A, A, • • • , A, • • •)] 

i=i 

= Pf(A5 Aj ■ ■ ■) + Pf(d-i, A) Aj di, • • •) • 

An algebraie identity of pfaffian^ 

Pf (d_i, A, A, A, • • • )Pf (■ ■ ■) = Pf (d-i, A, • • • )Pf(A, A, • • •) 

-Pf (d_i, A, • • • )Pf(A, A, • • ■) + Pf (d-i, A, • • • )Pf(A, A, • • •), 

implies 

(AAdi - 9i) X f = dsf X Afi'i - AA/ X S'* + Afi-i X dyf , 

whieh is aetually the first bilinear equation. 


dyds dy 
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The second bilinear equation is proved in a similar way as the one used by Iwao and Hirota^ 
in connection with a different system. We start from the r.h.s of the bilinear equation. 

l<fi<u<n 

. 27V 

= 4 E c,,J2^-iy^m{/3,,b,)Pf{dor--Mr--)mPu,b,)Pi{dor--,hr 

l</i,z/<n i,j 

27V . 

= ^c,,pm, k)pm, b,)Pi{do, • • •, ■ )pf(do, 

ij 

2A^ 

= (Pi" - Pf) )Pf(4, ■ • • , i;, ■ ■ • ) 

ij 

( 20 ) 

Next, the expansion of the vanishing pfaffian Pf( 6 j, do, • • ■) oii h yields 

2N 

J2(-iy+m{bi, 6 ,)Pf(do, ■■■ ,b„---) = Pf(do, • • •, d„ • • •), 
i=i 

which subsequently leads to 

2N 

5^(-l)*+^p-2pf(6,, 6 ,)Pf(do, • • • , • • • )Pf(do, •••,&„•••) 

2N 

= ^p“^Pf(do, • • • , di, • • • )Pf(do, • • • , 7i, • • •) • (21) 

i 

Similarly, we can show 

2N 

- ■ ■ ■ . ii. ■ • • )Pt(<io. • • ■ . ■ • •) 

id 

2N 

= ^p-2pf(do, • • • , d,, • • ■ )Pf(do, •••,&„•••)• (22) 

j 

Substituting Eqs. (l 2 TI)-(l 2 T]) into Eq. (l 20 l) . we arrive at 

^ 2N 

2 X] = ^a.^...)Pf(c^Q^... ^ 5 .^...). (23) 

l</i<z/<n i 


8 


Now we work on the l.h.s. of the second bilinear equation 

^ df Of 
ds^ ^ ds ds 

= Pi{d. 2 , do,--- )Pf(4, do,---)- Pf (d-i, do,--- )Pi{d.i, do,---) 

2N 2N 

= 5^(-l)Tf(d_2, a,)Pf{do, --- ,d„---) 5^(-l)^Pf(do, a,)Pf{do, --- ,dj,---) 

i=l j=l 

2N 2N 

- a,)Pf(do, --- ,d^,---) 5^(-iyPf(d_i, a,)Pf(do, --- ,d„---) 

i=l j=l 

2N 

= H W-2, a^)Pf{do, a,) - Pf(d_i, a,)Pf (ci_i, a,)] 

*,i=i 

xPf(do, --- ,di,--- )Pf(cio, --- ,dj,---) 

2N 

= [vi^ +Vi^V'J^] Pf(ai,aj)Pf(do,--- ,di, - --)Pf{do, - - - ,dj,---) 

The summation over the second term within the bracket vanishes due to the fact that 

2N 

Y^[-iy+3+^p-^p-T^Pi[ai,aj)Pi{do, --- ,di,- • •)Pf('^o, --- ,dj,---) 
ij=l 
2N 

= Y^{-iy^"^^p~^p~^Pi{aj,ai)Pi{do, •••,%,• • ■)Pf('^o, --- ,di,---) 

j.*=i 

2N 

= - Y^{-lY""^""^p~^p~^Pi{ai,aj)Pi{do, --- ,di,- • ■)Pf(c^o, •••,%,•••)• 

ij=l 

Therefore, 


dfdf 


2N 


= E (-1)*+^'^Vr^Pf(«*, a,)Pf(do, • • • , a„ • • ■ )Pf(do, ---,d„---) 


*j=i 


27V 


^(-l)*+Vi^Pf(c?o,--- ,di,--- 


2=1 


2N 


5^(-iyPf(a„a,)Pf(do,--- 

.i=i 


Further, we note that the following identity can be substituted into the term within bracket 


27V 


5^(-l)^Pf(a„a,)Pf(do,---,%,---) 

i=i 

= Pf(do, a,)Pf(- ■ ■) + (-l)*+'Pf (rfo, ---,k---) 


which is obtained from the expansion of the following vanishing pfaffian Pf(ai, do, • ‘ 
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Consequently, we have 

ds ds 

2N 


2=1 


5^(-l)‘+Vr'Pf (c?o, • • • , • • ■) [Pf (do, a,)Pf (• • •) + (-l)*+'Pf(do, ■■■Mr--) 

2N 

Pf (■ ■ ■ )Pf (d_ 2 , dor--) + Yl Pr'Pf (c^o, ---Mi,--- )mdo, ---,k---), 


2=1 


M f 

J A / y ^iiu9fi9u, 


5s 2 


(24) 


whieh ean be rewritten as 


, dfdf I ^ 

= 7T zZ <^i^^9iA9u- 


ds"^ 


ds ds 2 


(25) 


l<li<v<n 

The above equation is nothing but the seeond bilinear equation. Thus, the proof is eomplete. □ 


III. SEMI-DISCRETE ANALOGUE OF THE MULTI-COMPONENT SHORT PULSE 
EQUATION 


In this seetion, we attempt to eonstruet an integrable semi-diseretization of the MCSP equation 
dH). Firstly, we propose a semi-diserete analogue of bilinear equations dV]) 

1 


Ds{9ki ■ fk - 9k’ - fk+i) = dfc+i/fc + 9k’fk+i, j = 1,2,--- ,n, 


Sj) 


Sj) 


(i). 


Dlfk - fk = \ Mgtkk 


l< 2 <j'<ri. 

By introdueing a dependent variable transformation 

(h 

W _ i — 1 2 

~ fk' ’ ’ 

and a diserete version of the hodograph transformation 


(26) 


(27) 


Xk = 2ka-2{\nfk)s, t = s, 


(28) 


the seeond bilinear equation in (1261) is rewritten as 


(In/fc),, = - 


(*) (i) 
9M9k’ 


4 p 

l<i<j<n k 


1 

4 


E(d O') 

l<i<j<n 


(29) 
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From the discrete hodograph transformation, we can define an nonuniform mesh 

fk+l 


5k = Xk+i - Xfc = 2a - 2 In 


fk 


( 30 ) 


it then immediately follows 


dSk _ 1 


(31) 


., 0 ) 

l<2<ji<n 

from Eq. (|29|) . Next, dividing both sides by fk+ifk, the first bilinear equation in (|2^ can be 
calculated out by 

9kli,s 9k!s\ 9klifk,s - 9k^fk+i,s ^ ^ ^ 


fk+l fk 


fk+lfk 


fk+l fk I 


or 


a® 

/ 9k+i 9k 


^ , I 4k + _ f^\ 4k + ^ 

/fc+i fk / V A+1 fk ) \ fk+l fk 


\ fk+l fk 

which is recast into 

\ fk+l fk I \ \ fk 

With the use of Eqs. (1271) and (l30l) . we finally arrive at 

c((4+i - uf) 1 


= I a- + ^ 


\ fk+l fk 


(32) 


= ^(a^fc+i - (4+1 + 4*^)- 


ds 


(33) 


Eqs. (|3T]) . (1331) constitute the semi-discrete analogue of the MCSP equation. We summarize the 
results by the following Theorem. 

Theorem 2. The bilinear equations ( 1231) yield a semi-discrete multi-component short pulse equa¬ 
tion 


< 


^ d{u^k+i ~ “k) 1 r / (i) (jT 


ds 

dSk _ 1 

ds 




(0 (i) (*) (i) 

c*j i4;i4+i 4 


(34) 


through dependent variable transformation 


(h 

~ * = 3) 2, • • • , n , 

fk 


and discrete hodograph transformation 


Xk = 2ka-2{\nfk)s, t = s, 


where 5k = Xk+i - Xk- 
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To assure its integrability, we provide its multi-soliton solution in terms of pfaffians by the 
following theorem. The elements of the pfaffians are defined as follows: 


ai.afk - , Vi \k)T\ \k) , Fi{aiffij)k - bij , 

(35) 

Pf(6„ bfk = 7 2 > (b, e B^, b, e B,) , 

4: Pi -pj 

(36) 

Pf{di,afik = Tf\k), Pf{a„d^)k = ^‘f\k + 1), 

(37) 

f 1 bj e Ba 

Pi{b,,fi,)k = { ' ^ 

[ 0 bj ^ Bfj, 

(38) 

Pf{d,,d^)k = l, Pf{d_,,d% = -a, 

(39) 


where 


I ^ \ t_-1. 


pr{k) = p^ 


ii=Pi s + 6 


iO 


I - api 

Here i,j = 1, 2, • • • , 2N, - ,n and k,l are arbitrary integers. Other pfaffian elements 

not mentioned above are all zeros. Note that (p\^\k) has the following property 

v,r\k+i)-^P{k) ^ ^ _ 

a 

which is used in the proof of the theorem. 

Theorem 3. The bilinear equations ([2^ have the following pfaffian solution 

fk = Pf(ai, • • • , a2N, bi, - ■ ■ , &2Af)fc , (40) 

= Fi{do,flj,ai, ■ • • ,a2Ar,fei, • • • , 62 iv)fc , (41) 

where j = 1, 2, • • • ,n, the elements of pfaffians are defined in eqs. 

Proof Since 

ds 


Fi{ai,aj)k = ^f\k)(p) ^\k) - ip\ ^\k)ipf\k) = Pf(ci_i, do, a*, , 


Pf(ai,aj)fc+i = Pf(ai,aj)A: + T?\k + l)pf\k) - pf\k)pf\k + 1) 

= Pf(do,d*^,ai,aj)fc, 
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= vf\k + + 1) - vt^\k + l)ipf\k + 1) 

-a (^Pf(ai, aj)k + ipf\k + l)ipf\k) - ipf\k)ipf\k + 1)) 
= -aPf(ai, aj)k + (pf\k + 1) (^ip^~^\k + 1) - aipf\k)^ 
-ipf\k + 1) {ip^~^\k + 1) - aipf\k)^ 

= -aPf{ai, aj)k + ^f\k + 1 ) (ai^f\k + 1 ) + (p^~^\k)^ 
-ipf\k + 1) (aipf\k + 1) + 

= -aPi{ai,aj)k + <pf\k + l)(pY^\k) - (p^~^\k)(pf\k + 1) 
Pf((i_i, d , flj, cij^k ) 


we have 


dJk = Pi{d_udo,---)k, 


fk+i = Piido,d^---)k, 


{ds - a) fk+i = Pf(ci-i, d^ ■■■)k. 


Furthermore, we can verify 


dsgi^^ = d, 



-l)*Pf(do,a.)fcPf(A 


^l1 


2N 

= (('9sPf(cio, a,)k) PfW^, ■■■ ,d„---)k + Pf{do, a,)kdsPf{/3^, • • • , a,, • • • )fc) 

i=l 

2N 

= ^(-1)* (Pf(ci-1, ai)kPi{/3^„ ■■■ ,di,---)k + Pi{do, ai)fcPf(/3^, d-i, do,-" , a*, • • • )k) 

i=l 


= Pf{d-i,/3^, ■■■)k + Pi{do,/3^,d-i,do, ■■■)k 

= Pf{d,,/3„---)k, 
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2N 

= 5^(-l)Tf (rfo, a,)k+iPiW„ • • • , a,, • • • )fc+i 

i=l 

27V 

= 5^(-l)*-lpf(ci^ a,)fcPf (/9^, do, d\---,di,---)k 

i=l 

27V / 2-1 

= ^(-i)-‘pf(o!‘,OiK Pt(/3„■ • ■,a,.■ ■ ■ )t + ■ “J"■ ■ ■ ■ 

i=i V j=i 

2N 

j=i+l 

= Pf(o!‘,/3„,-'-)i. 



2N 

(a. - a)gi% = ^(-1)-' ((a. - <i)Pf(a‘, oOO Pt(4. • • ■. a,. ■ ■ ■ 

2 = 1 

+ (-1)*“^^"^ (^aPf(ai,aj)fc) Pf(/g^,rf'', • • • ,ai, • • • ,aj, • • •)/= 

l<i<j<2Af 

2N 

= {Pf{d-1, a,), + aPf(rfo, a,)k) PiW„ • • • , a„ • • • )fc 

2=1 

+ c?o, a*, aj)fcPf (/?/,, d^,--- ,di,--- ,%,••• )fc 

l<i<jr<2Af 

= Pi{d.i,l3^,do,d\---)k. 

Therefore, an algebraie identity of pfaffian 

Pf (d_i, 13^, do, d\--- )fcPf (• ■■)k = Pf(d-i, do,--- )kPi{dy (3^,--- )k 
-Pf (d_i, • • • )kPi{do, /3„---)k + Pi{d-i,/3^, - - - )fcPf (c?o, dy---)k, 

together with above pfaffian relations gives 

(ds - a)giy^ X fk = gi+i X dsfk - {ds - a)fk+i x 4+i + dsg^^ x fk+i, 

whieh is nothing but the first bilinear equation. The seeond bilinear equation ean be proved in a 
similar way as in the continuous case. □ 

Remark 4. Bilinear equations can be viewed as a Backlund transformation of bilinear equa¬ 
tions (IZl), which yield the MCSP equation. In other words, if fk and g^^ satisfy (0), so do fk+i and 
g^l^^ . Based on ^2M . we propose a semi-discrete analogue of the MCSP equation. The integrabil- 
ity of the semi-discrete MCSP equation is guaranteed by the existence of N-soliton solution. 
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Finally, let us show that in the eontinuous limit, a —)■ 0 (5^ —)■ 0), the proposed semi-discrete 
multi-component short pulse equation recovers the continuous one (flSl) . The dependent variable 
u is regarded as a function of x and t, where x is the space coordinate of the k-th lattice point and 
t is the time, defined by 

k-l 

X = xq + 5j, t = s, 
j=0 


where 6j = Xj+i — Xj. In the continuous limit, a ^ 0 (5^ —>■ 0), we have 


-[u 


(«i+i + “fcO ^ 




dju 


(d 

/c+1 


U 




-)• Ui 


dx dxo ^ dSj 
ds ds ^ ds 

j=0 


dx 


k-l 


0 


ds 


hence 


2 E 

l<fj.<iy<n j=0 

l<^<v<n 


\l<lJ,<y 


ip) i'^) ik) M 


I 9x 5 


dXn 


Thus 


where the origin of space coordinate Xq is taken so that cancels ^ J2i</_i<u<n “o 

the first semi-discrete multi-component SP equation converges to 

1 


d^ dt 


( ^ ^ I dx I Ui Ui , 

\l<fi<u<n / / 


which is exactly the MCSP equation (jH). 


IV. TWO-COMPONENT SHORT PULSE EQUATION 

Since the two-component short pulse equation is of particular importance for applications in 
nonlinear optics, we provide a detailed study for this two-component system, together with its 
semi-discrete analogue in this section. For the continuous case of n = 2, we can take Ci 2 = 1 
without loss of generality and arrive at the following two-component system— 

Uxt = U+]^ {uVUx)x , ( 42 ) 

Vxt=V + ]^ {uVVx)x , ( 43 ) 
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where u = ui,v = U 2 . Furthermore, if we assume m is a eomplex-values function and impose a 
complex conjugate condition v = u, where u means the complex conjugate of u. Eq. (|4^ leads to 
a complex short pulse equation studied in22^ 


^xt 





(44) 


Since the complex short pulse equation (l44l) is a special case of two-component system (I42l)-(l43l). 
its A^-soliton solution can be obtained from the iV-soliton solution of system (l42l) - (l4^ by requiring 
/ = /, These requirements can be achieved by putting = pj and ,^^0 = where 

k= j + N,j = 1,2,■■■ ,N. 

In particular, the tau-functions for one-soliton solution (N = 1) are found to be 


1 {pipif 

^ 4(pi+pi)2 


9 = 


(45) 

(46) 


Let Pi = piR + ipij, and we assume pm > 0 without loss of generality, then the one-soliton 
solution can be expressed in the following parametric form 


u = ?^e'^i^sech {gm + pio) , 


|pil 2 


x = y- 


"^VlR 

\P1? 


where 


(tanh {giR + r/io) + 1) , t = s. 


, PiR Pii 1 \Pi\ 

ViR = PirV + 1 pu = PiiV - , pio = In -- 

\Pl\ Pi 4piH 


(47) 

(48) 

(49) 


Eq. (l47l) represents an envelope soliton of amplitude 2piR/\pi\‘^ and phase pu. The details 
analysis concerning its property was carried out in^i. In summarythree types can be classified. 

• when \pir\ < \pii\, it is a smooth soliton solution, which is similar to the envelope soliton 
solution for the nonlinear Schrddinger equation. 

• when \pir\ > \pii\, it is a loop soliton solution, which admits multi-valued property. 


when \pir\ = \pii\, it is a cuspon soliton, this is a case which divides the single-valued and 
multi-valued solution. 
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A. Semi-discrete two-component system 


Based on the results in previous section, we have an integrable semi-discrete analogue of two- 
component system (I42l)-(l4^ 


Ufg'j ^ r / \ 

“f '^kjj 

dSj^ 1 , , 

-TT = --[Uk+lVk+1 -UkVk), 

at 2 


which admit the following iV-soliton solution in parametric form 


(50) 


( 1 ) ( 2 ) 
Jk Jk 


(51) 


and hodograph transformation 


Xk = 2ka-2{\nfk)s, t = s, 


(52) 


with 


fk — Pf(ai5 • • • , o,2N, bir '' ) b^, Ci, • • • , cisf)k , 


Qk P (*(/^^5 ^ 0 ? ^ 1 : * * * : ^2N 'I b\i ‘ * * 5 bp^^ ‘ ‘ * 7 ^N^k } ^ f 7 2 


(53) 

(54) 


where the elements of pfaffians are (others not mentioned are all zeros) 


Pf(aj,aj)fc = ipf\k)ipf\k), Pf{ai,bj)k = Sij, 

Pi + Pj 


Pi{h,c,) 


1 {PiPN+jf 
^pf-PN+j ’ 


Pf((2j, Cj')k , 


Pf(6„ A) = Pf(Q,/?2) = 1, Pf(do,a.)fc = V^f\k). 


By imposing complex conjugate conditions pk = Pj, ^ko = ^jo (k = j + N, j = 1, 2, • • • , N), it 
then follows fk = fk, g^k'^ = g^k'^, thus Uk = Vk, which leads to a semi-discrete analogue of the 
complex short equation (l44l) 


{ d{Uk+l Uk) _ Ir / . \ 

'^Ok\Uk+l + Uk) , 


( 55 ) 
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Its iV-soliton solution immediately follows from the iV-soliton solution of Eqs. (l5^ - (l54l) under 
complex conjugate conditions mentioned above. In what follows, we list the one- and two-soliton 
solutions. 

One-soliton solution: The tau-functions for one-soliton solution to Eq. (1551) are 

/fc = -1 - (k), 9^^ = (56) 


with 



f 1 + api'\’' p-h+$io (0)n ^ 




where pi = p 2 , ^lo = ^ 2 o- Similar to the continuous case, if pi = pm + ipu, we then arrive at the 
one-soliton solution of semi-discrete complex short pulse equation (1551) 


Xk =2ka— 


where 


a UA I 

Ok = kdi + - — -s. 


= e'^'=sech(6'fc-f (lo), 

bil 

(57) 

1 12 (tanh(6'fc -f Oq) + 1) , 
bil 

(58) 

|pi|^ 1-opi 

(59) 


\pi? ’ 

In Eig. 1 (a)-(c), we illustrate the envelope soliton for pi = 1 + 1.5i, 1 -f i, 1 -f 0.5i, which 
correspond to the smooth, cuspon and loop solition, respectively in the continuous case. Two- 
soliton: The tau-functions for two-soliton solutions to semi-discrete complex short pulse equation 
(1551) can be obtained by 

fk = l + 013(^13 (fc) + aiiPlhk) + 023^^23^ (fc) + «249524 (^) 

-f 16012023014024 {Pl^ - P2^Y - PI^Y pfhk){k)pfl{k) , 


9^Y = YYik) + YYik) + 4 (^0i3023P® (fc) + aiia2ipf\k)^ {p^ ^ - P 2 Puik) ■ 

where 

Pi -Pj 

a-ij = —^-, Pi=P3, P2=Pi 

Pi + Pj 

and .^10 = ^30 and ,^20 = 'f 40 - In Eig. 2 (a)-(d), we show the process of interaction between a 
smooth envelop soliton and a cuspon envelop soliton with pi = 1 + 1.5i, p 2 = 1 + 1 , ^10 = —15, 
60 = -25. 
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FIG. 1. Envelope soliton for the semi-discrete complex short pulse equation (a) smooth soliton with pi = 
1 -|- 1.5i, (b) cuspon soliton with pi = 1 -|- i, (c) loop soliotn with pi = 1 -|- 0.5i. 

V. CONCLUDING REMARKS 

We have derived an integrable semi-diserete analogue of the multi-eomponent short pulse equa¬ 
tion proposed by Matsuno^i based on a Baeklund transform and Hirota’s bilinear method. We find 
its A^-soliton solution in terms of pfaffians and prove it. Moreover, a eomplex short pulse equa¬ 
tion, which possess smooth, cuspon or loop type envelop soliton, is proposed and its semi-discrete 
analogue is constructed as well. We conclude the present paper by the following remarks. 

• The A^-solution for multi-component short pulse equation given in the present paper agrees 
with the one given by Matsuno in^. This solution is a benchmark for the study of soliton 
interactions. 

• Similar to our previous results^^i^^, the semi-discrete multi-component short pulse equa- 
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X X 


FIG. 2. Two soliton solution for the semi-discrete complex short pulse equation; solid line: |tt|, dashed line: 
Reiu). (a) i = 0; (b) f = 45.0; (c) t = 60.0; (d) t = 110.0. 

tion proposed here can be served as an integrable numerical scheme, the so-called self- 
adaptive moving mesh method, for the numerical simulation of multi-component short pulse 
equation, as well as complex and coupled complex short pulse equation. However, it is ob¬ 
viously beyond the scope of the present paper, we would like to report our results in this 
aspect in a forthcoming paper. 


• The integrable fully discretization of the multi-component short pulse equation is a further 
topic deserve to study. 
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